Majorana Fermions Under Uniaxial Stress in Semiconductor-Superconductor 

Heterostructures 
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Spin-orbit coupled semiconductor nanowires with Zeeman splitting in proximity contact with 
bulk s-wave superconductivity have recently been proposed as a promising platform for realizing 
Majorana fermions. However, in this setup the chemical potential of the nanowire is generally 
pinned by the Fermi surface of the superconductor. This makes the tuning of the chemical potential 
by external electrical gates, a crucial requirement for unambiguous detection of Majorana fermions, 
very challenging in experiments. Here we show that tunable topological superconducting regime 
supporting Majorana fermions can be realized in semiconductor nanowires using uniaxial stress. For 
n-type nanowires the uniaxial stress tunes the effective chemical potential, while for p-type systems 
the effective pairing may also be modified by stress, thus significantly enhancing the topological 
minigap. We show that the required stress, of the order of 0.1%, is within current experimental 
reach using conventional piezo crystals. 

PACS numbers: 71.10.Pm, 03.67.Lx, 74.45.+C, 74.78.-w 



Majorana fermions (MFs) are quantum particles which 
are their own antiparticles^— . MFs are not only of fun- 
damental interest because of their non-Abelian exchange 
statistics, but also may serve as building blocks for fault- 
tolerant topological quantum computation 4 ' 5 . In the past 
few years, the possibility of realizing MFs using quasi- 
particles in exotic solid state systems^— has generated 
a lot of excitement in the condensed matter community. 
In particular, it has been proposed recently that MFs 
can be generated using a heterostructure consisting of 
two very conventional materials: an s-wave superconduc- 
tor and a semiconductor thin film/nanowire with strong 
spin-orbit coupling 1 1-20 . This proposal, following on the 
earlier idea proposed in cold atoms systems 2 ^, has at- 
tracted widespread theoretical interest as well as serious 
consideration in experiments. Very recently, following 
the theoretical proposals, some preliminary experimental 
signatures^— which may be related to the existence of 
MFs have been observed, although the unambiguous de- 
tection of MFs still remains an outstanding experimental 
challenge. 

In the proposed semiconductor heterostructures, MFs 
only exist for topological superconducting states within 
a certain constrained parameter regime 13 . Therefore, it 
is crucially important to have the ability to tune various 
physical parameters in experiments. Since the spin-orbit 
coupling strength and the size of the nanowire cannot 
be tuned after the samples are fabricated, and the Zee- 
man field V z is generally limited to a narrow window 
{\V Z \ ~ A, A is the superconducting pairing) because of 
its possible depairing effect, it is essential to be able to 
tune the chemical potential to the correct level for the 
realization of the topological state. For semiconductor 
nanowires, the chemical potential can be tuned by using 
suitably placed external electrical gates. Unfortunately, 
the same technique may not work well in the proposed 
Majorana system because the nanowire is in proximity 



(a) Model A 



Piezo Crystal 



(c) Pinned Chemical Potential by The 
Fermi Surface of Superconductor 
/ 
'// 



S-wave Superconductor 



(b) Model B 



Piezo Crystal 



S-wave Superconductor 




' Tuned 
by stress 



Band structure 
in nanowire 



B // z Fermi Sea of S-wave 
Superconductor 



FIG. 1: (Color online). Engineering Majorana fermions via 
uniaxial stress. The stretching direction is along x in model A 
(a) and z in model B (b) . The stress is assumed to be provided 
by conventional piezo crystals, (c) shows the pinning of the 
chemical potential in semiconductor by the Fermi surface of 
the s-wave superconductor. However, the total band structure 
of the semiconductor can still be shifted by uniaxial stress. 



contact with the superconductor which has an extremely 
high carrier density. The pinning of the chemical poten- 
tial of the nanowire by the Fermi surface of the super- 
conductor thus poses a major challenge in experiments^. 
In the recent Delft experiments 2 ^ it has been found that 
very large electric gate voltages (~ (10 4 — 10 5 )A) are 
required to tune the topological quantum phase transi- 
tion, which likely provides strong evidence for the pin- 
ning of the nanowire chemical potential. The main mo- 
tivation of the present work is to provide an alternative 
method to overcome this experimental difficulty which 
may greatly facilitate future experimental searches of 
MFs in semiconductor-based heterostructures. 

In this Rapid Communication we show that the topo- 
logical superconducting regime in both n- and p-type 
nanowires can be externally tuned using uniaxial stress, 
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which can be generated and controlled even by conven- 
tional piezo crystals. The uniaxial stress can modify the 
band structure of the nanowires slightly after the het- 
erostructure system is fabricated, remarkably leading to 
a topological transition from the trivial to the topolog- 
ical superconducting state with MFs at the wire-ends. 
With the experimentally accessible strength of the uni- 
axial stress, the effective chemical potential can be tuned 
about 42 meV for electrons, and 5-20 meV for hole 
levels. Moreover, for the p-type systems, the uniaxial 
stress may also significantly enhance the minimum topo- 
logical energy gap (minigap) that protects the MFs from 
thermal excitations. The newly added elements for gen- 
erating the uniaxial stress can be effectively integrated 
into the design of semiconductor devices using modern 
nanotechnology. Therefore, our proposed scheme can go 
a long way in facilitating the realization and detection of 
Majorana fermions in semiconductor quantum wire het- 
erostructures and the eventual implementation of topo- 
logical quantum computation. 

Our basic setup for experiments is illustrated in Fig. [1]} 
and[T|3. The semiconductor nanowire (e.g., InSb, InAs, 
etc.) is in proximity contact with an s-wave supercon- 
ductor. The uniaxial stress applied on the semiconduc- 
tor nanowires can be generated using nano-ferroelectric 
materials or by simply gluing the nanowire tightly to the 
surface of piezoelectric crystals^ 9 -"— such as the piezo- 
electric lead zirconic titanate (PZT) ceramic stack. The 
stretching direction of the piezo crystal can be chosen ei- 
ther along x (model A) or along z direction (model B). 
The strain tensor can be determined as, 

F (o) = _ F F {a) = Jo) = 2C i2 

On 



where e = (1 — a/ao) defines the relative changes of 
the lattice constant along the corresponding crystallo- 
graphic directions. Here ao and a are the equilibrium 
and distorted lattice constants, respectively, e > 0(< 0) 
corresponds to compressive (tensile) stress. In experi- 
ments the sign of e can be controlled by the voltage bias 
across the piezo crystals^— . Cn and C\i are the elas- 
tic stiffness tensors. The superscripts in Eq. (fTJ) repre- 
sent the two different models shown in Fig. [1] Note that 
e ~ V/y, where V is the stress and y is Young's modu- 
lus. Using typical values for y ~ 100 GPa and V ~ 100 
MPa, we see that e ~ 0.1%. Such a small strain can be 
provided using conventional piezo crystals. For model A, 
e up to 0.11% has already been realized in experiments, 
and in principle, e up to 0.6% can be achieve d 30 ' 31 . For 
model B, it is more suitable to provide compressive stress 
along the z direction, and there is no limitation on the 
maximum e because the compressive stress is not lim- 
ited by the gluing technique. We assume |e| < 0.3% 
the most probable regime that can be accessed in exper- 
iments. Since the lattice deformation is very small, the 
uniaxial stress has negligible effects on the superconduc- 
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FIG. 2: (Color online) (a) Band structure of InSb nanowires. 
The solid horizontal dash-dotted line represents the possible 
pinned chemical potential when placed in proximity to the 
superconductor, (b) The chemical potential as a function of 
carrier density, a ~ 0.2 eV-A, m = 0.013mo, Ao = 0.5 meV, 
Vz = 1.0 meV. Parameters are from Ref. |3y. 

tor transition temperature as well as the s-wave pairing 
symmetr y 34 ' 35 , two properties that are crucial for the 
generation of Majorana fermions. 

The effective Hamiltonian for n-type nanowire under 
uniaxial stress reads a o 38 ' 39 , 

k 2 

H = [— + a(p x a) z + V z a z + a c Tr(e)] - fi F . (2) 
2m 

where a c denotes the deformation potential of the con- 
duction band, a is the spin-orbit coupling strength, and 
V z is the external Zeeman field induced by the magnetic 
field, hf is the true chemical potential of the semicon- 
ductor that is pinned to the Fermi surface of the super- 
conductor. The uniaxial stress does not change the band 
structure, but shifts the effective chemical potential to 
/i = /if — a c Tr(e), through which the topological region 
V 2 > A 2 + jjpl&lS. may be achieved. In Fig. [5h., we plot 
the typical band structure of free electrons in nanowires 
in a single transverse confinement band and in Fig. [2Jd 
we plot the corresponding chemical potential as a func- 
tion of carrier density. The topological superconductivity 
with MFs can be achieved when the chemical potential 
falls in the small window in Fig. ^Bp : in which case the 
system cuts only one Fermi surface. The small window in 
the parameter space greatly limits the flexibility for the 
experimental observation of MFs. For InSb, a c ~ —6.94 
eV, and |e| < 0.3%, we estimate a c Tr(e) ~ ±21 meV, 
the same magnitude as the energy difference between the 
chemical potential of the nanowire and the conduction 
band minima in the Delft experimental. We see from 
Fig. [2b that such a huge change of the effective chemical 
potential can change the carrier density by about 1-2 
orders of magnitude. Thus, for a wide range of carrier 
density, the nanowire can always be tuned to the topo- 
logical regime in the experiment. 

The Hamiltonian for p-type nanowires (assuming the 
sizes along the y and z directions are L y and L z , respec- 
tively) under uniaxial stress reads as^i, 

H = Hkl + Hbp — [i-Fi (3) 
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FIG. 3: (Color online), (a) Typical band dispersions of the 
heavy hole bands for different values of strain e in model A. 
(b) Energy of heavy hole bands at k — as a function of 
strain e for models A and B. Other parameters are 71 = 34.8, 
72 = 15.5, a = 2.0 eV-A, V z = 1.5 meV, L z = 14 nm, and 
L y = 10 nm, fi = 0. Parameters are from Ref. |36| . 



count. Because the uniaxial stress only shifts the effec- 
tive chemical potential for the n-type of semiconductor, 
we consider only p-type nanowires. In the nanowire het- 
erostructures, the superconducting order parameter can 
be induced to the nanowire through proximity effec t 42 ' 43 , 
yielding H sc = X)m=i,| / drA m^ m ^- m - The corre- 
sponding BdG Hamiltonian can be written asi 9 - 



H 



BdG 



iJlD A4 

A 4 -7 f ^iD7 



(5) 



in the Nambu spinor basis ^ = (ij}, r yip') T . Here H 
J dydzi/}*ip*Hip z ip y , 7 = i(I ® a x )r y 



I and t v are 



Pauli operators, and A 4 =diag(As, Ai, Ai, A3). 

222 2 

The topological parameter regime for MFs can be ob- 
tained by the topological inde x 44 ' 45 



where F KL = ^+ 5 ^ V 2 - 72 (V-J) 2 - M (JxV) z + 14Jz 
is the Kohn-Luttinger Hamiltonian, J is the total angular 
momentum operator for the spin-3/2 holes and 71 and 72 
are Luttinger parameters. The second term describes the 
Bir-Pikus model^ii 



M = sign(Pf(r(0)) • Pf(F(-))) 
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e yy ) : with a v and b the deformation 
potentials of the valence bands. Notice that P e , Q e and 
R £ have totally different roles to the band structures of 
holes in nanowires. P £ shifts the global band structure, 
while Q e increases or decreases the splitting between the 
heavy hole (HH) and light hole (LH) bands. The non- 
zero R e may greatly enhance or suppress the coupling 
between HH and LH, thus modifying the effective pair- 
ing strength. In contrast to the case of electrons, the two 
models A and B yield totally different results. 

The modification of the band structure of p-type 
nanowires due to uniaxial stress is shown in Fig. 3a for 
different values of e. Here we only plot the two HH bands 
because the LH bands are separated by a large energy gap 
(~ 100 meV) induced by the confinement. At k = 0, two 
HH bands are split by a small Zeeman field Vz- When 
the chemical potential lies in the Zeeman gap, the system 
has only a single Fermi surface, yielding topological su- 
perconductivity and the associated MFs under suitable 
conditions. We see from Fig. 3 that by tuning the uni- 
axial stress, we can shift the bands of the semiconductor 
up or down so that the Fermi level of the superconductor 
can lie in the Zeeman gap, yielding topological supercon- 
ductivity. In Fig. 3b, we plot Eq, the energy of the HH 
bands at k = 0, against e. Within the experimentally 
accessible regime, Eg can be tuned by ±3 meV for model 
B and ±10 meV for model A. 

To obtain concrete parameter regions of the MFs, the 
superconducting order parameter need be taken into ac- 



where PffT) refers to the Pfaffian of the matrix T — 
—iHBdG(k){T y ® 7), a is the lattice constant. Ai = +1 
(-1) corresponds to the topologically trivial (nontrivial) 
superconducting states without (with) MFs. Note that 
for sufficient large k all the eigenvalues of F are dom- 
inated by the k 2 terms, yielding sign(Pf(r(^))) = 1. 
The Pfaffian at k = can be derived analytically, yield- 
ing M = sgnJ", with T = f - fiV z 2 + ^V z 4 , f = 
( / 2 2 + A3A.-/3 1 2 -/3 2 2 ) 2 + ((A | -A.)/2+/3 1 (A | + A.)) 2 , 
/1 = (10/i 2 + 10/3? + 16^/2 + 9A| + A| - 6/3f)/4, 

Pi = 7r 2 7 2 (L- 2 -L- 2 /2) + Q e , /? 2 = V37r 2 72 L- 2 /2 + i? e , 
and p, = fi + -f^lly 2 + L- 2 /2) 4e . The boundary for 
topological phase transition is determined by J 7 = 0. 
Generally, the magnitudes of A 3 and Ai are not es- 
sential for the topological quantum phase transition (but 
the relative sign is important) 19 . Henceforth we only con- 
sider two different possible cases (I) A3 = Ai = A and 
(II) A3 = — Ai = A. For other values of A3 and Ai, 

2 2 t 2 2 

the results are similar. 

In Fig. 2k and Hp, we plot the boundary between the 
topological and non-topological superconducting states. 
In these figures, we assume that without uniaxial stress 
the chemical potential lies in a regime which requires a 
large Zeeman field for realizing the topological supercon- 
ducting state. By applying the uniaxial stress the re- 
quired critical Zeeman field can be significantly reduced. 
In case (I) in Fig. the critical Zeeman field can 

even approaches zero for model A. We have also veri- 
fied that for a wide range of parameters (/z, L Zl L y • • ■ ) 
similar features can always be found. For case (II) in 
Fig. Hp the required Zeeman field can be reduced to 
around 1 meV (B z = 0.3 T for g\ - 50). Generally 
for case (II), the minimum required V z ~ pA, where 

p = 2y/T+pJ7M/& + V^+W/M) e [2/3, 2). We see 
that there are also a wide range of parameters that en- 
able us to achieve the minimum required Zeeman field 
at p — 2/3 via uniaxial stress. To further verify that 
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FIG. 4: (Color online). The parameter regimes for the exis- 
tence of Majorana fermions. (a), (b) show the results for case 
(I) and case (II), respectively, see text for details. The plus 
sign (+) corresponds to V z without uniaxial stress. For case 
(I) we use L y ~ 10.0 nm, L z ~ 14.0 nm, for case (II) we use 
L y ~ 14.2 nm, L z ~ 9.7 nm. The chemical potential in all 
the figures are pinned at \i = —96.0 meV. In (b) the shaded 
regime corresponds to the lower bound of V z /A = 2/3. The 
corresponding mini-gap and lowest non-negative energy level 
obtained by numerically solving the BdG Eq. ([5]) are pre- 
sented in (c) and (d), respectively, for V z = 1.5 meV. The 
other parameters are the same as Fig. [3] 



the right regime of each curve in Fig. 0^ and0}) are in- 
deed the topological superconducting regime we plot the 
mini-gap (solid line) and the lowest non-negative energy 
level (dashed line) as a function of uniaxial stress e in 
Figs. I3J; and 0)1, respectively. In the topological super- 
conducting state, the zero energy MFs indeed exist with 
large minigaps around several Kelvin. We have also con- 
firmed that the corresponding wavefunctions of MFs are 
well localized at the two ends of the nanowirc. 

The difference between the two types of uniaxial stress 
in models A and B can be understood by projecting the 
Hamiltonian to the lowest two HH bands ^ 41 ' 46 ! 47 , which 
yields the effective pairing at k —> 0, 



A eff = (As -kAi)/(k- 

x 2 2 ' ' v 

WWi/M 2 + 1 - Nhf 



i 



(7) 



where k = + 1 - fa/faf € (0,oo). When 

e — 0, k only depends on the size of the nanowire, thus 
cannot be tuned. However, when the uniaxial stress is 
applied, k can be tuned in a considerably wide range. 
For model A, the off-diagonal term R e ^ 0, thus /32 can 



approach zero with a properly chosen strain £. For the 
parameters used in Fig. 2] we find that the effective pair- 
ing increases (decreases) monotonically as a function of e 
for model A (B), thus for model B, we observe significant 
enhancement of the mini-gap (~ 30%) in Fig. [4ji. The 
maximum increase of the mini-gap can be obtained by 
optimizing different physical parameters. 

Finally, several remarks are in order. First, the same 
idea discussed above for a single band model can be 
straightforwardly extended to the multiband case. Us- 
ing the diameter of nanowire from Refill, we estimate 
the band spacing for electron (hole) is ~ 6 (~ 3) meV, 
with ~ 4 bands occupied. Thus we expect that the stress 
can tunc the effective chemical potential of both elec- 
trons and holes to the topological regime even though 
the initial value corresponds to an even number of Fermi 
surfaces. Second, our proposal here can also be used to 
engineer MFs in the vortex core of semiconductor quan- 
tum wells. For electrons the tuning of the band struc- 
tures is exactly the same as in Eq. [2j For holes there are 
some qualitative difference from Fig. 3 since the confine- 
ment along the y direction is relaxed. As a consequence, 
/?2 = when £ = 0, thus R £ play a more significant role 
in the determination of the minigap of MFs. Third, we 
have also checked the validity of our proposal for InAs 
nanowires, and similar features have been found. How- 
ever, for InAs, we note that££ a c ~ —5.17 and b ~ —1.00, 
which are smaller than their counterparts in InSb, thus 
a slightly larger stress is required. 

To conclude, due to the proximity effect between 
nanowires and a bulk superconductor, the chemical po- 
tential of the nanowire is generally pinned by the Fermi 
surface of the superconductor. Consequently, tuning 
the chemical potential of nanowires via electrical gates 
to bring it in the topological regime is inefficient in 
this setup. We show that this crucial obstacle can be 
overcome using experimentally accessible uniaxial stress 
which modifies the band structure slightly, leading, re- 
markably, to a transition from non-topological to topo- 
logical states with MFs. The newly added elements for 
generating uniaxial stress can be effectively integrated 
into the design of semiconductor devices using modern 
nanotechnology. Therefore our scheme can be used for 
the realization of topological Majorana fermion excita- 
tions in semiconductors and the implementation of topo- 
logical quantum computation. 
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